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ABSTRACT 



We present the first results of AKARI Infrared Camera near-infrared spec- 
troscopic survey of the Large Magellanic Cloud (LMC). We detected absorption 
features of the H2O ice 3.05/im and the CO2 ice 4.27/im stretching mode toward 
seven massive young stellar objects (YSOs). These samples are for the first time 
spectroscopically confirmed to be YSOs. We used a curve-of-growth method 
to evaluate the column densities of the ices and derived the CO2/H2O ratio to 
be 0.45±0.17. This is clearly higher than that seen in Galactic massive YSOs 
(0.17±0.03). We suggest that the strong ultraviolet radiation field and/or the 
high dust temperature in the LMC may be responsible for the observed high CO2 
ice abundance. 

Subject headings: circumstellar matter — stars: pre-main sequence — ISM: 
abundances — dust, extinction — ISM: molecules — Magellanic Clouds 



1. Introduction 



Properties of extragalactic young stellar objects (YSOs) provide us important informa- 
tion on the understanding of the diversity of YSOs in different galactic environments . The 
Large Magellanic Cloud (LMC), the nearest irregular galaxy to our Galaxy (~50kpc; lAlves 
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20041 ), offers an ideal en vironment for this study since it holds a unique metal-poor environ- 
ment (iLuck et al.lll998l ). Because of its proximity and ne arly face-on geometry, various types 



of su r veys have been p erformed toward the LMC (e.g., 
20061 : iKato et al.l 120071 . and references therein) . 



Zaritsky et al.ll2004j : iMeixner et al. 



An infrared spectrum of YSOs shows absorption features of various ices which are 
thought to be an important reservoir of heavy elements and complex molecul es in a cold 



environment such as a dense molecular cloud or an envelope of a YSO (e.g., IChiar et al. 



19981 : iNummelin et al.lboOll : Iwhittet et aPboO?! : ISooeert et aPboOsh . These ices are thought 
to be taken into planets and comets as a result of subsequent planetary formation activity 
(lEhrenfreund fc Schuttd l2000l ). Studying the compositions of ices as functions of physical 
environments is crucial to understand the chemical evolution in circumstellar environments 
of YSOs and is a key topic of as trophysics. H2O and CO2 ic e s are ubiquitous and are majo r 
components of interstellar ices (Ivan Dishoeck fc Blakdll998l : iBoogert fc Ehrenfreundl l2004j ) . 
Since the absorption profile of the ices is sensitive to a chemical composition of icy grain 
mantles and a thermal history of local environments, the ices are important tracers to inves- 
tigate the properties of YSOs. However, our knowledge about the ices around extragalactic 
YSOs is hmited because few observations have been performed toward extragalactic YSOs. 
Therefore infrared spectroscopic observations toward YSOs in the LMC are important if we 
are to improve our understanding of the influence of galactic environments on the properties 
of YSOs and ices. 

AKARI is the first Japanese sate llite dedicated to an infrared astronomy launched in 
February 2006 (IMurakami et al.l 120071 ). We have performed a near infrared spectroscopic 
surve y of the LMC usin g a powerful spectroscopic survey capability of Infrared Camera 
(IRC; lOnaka et al.l 120071 ) on board AKARI. In this letter, we present 2.5-5;um spectra of 
newly confirmed YSOs in the LMC with our survey, and discuss the abundance of H2O and 
CO2 ice. 



2. Observations and Data Reduction 



The o bservations rep orted here were obtained as a part of the AKARI IRC survey of 
the LMC (jita et al.ll2008l ). An unbiased slitless prism spectroscopic survey of the LMC has 
been performed since May 2006. In this survey, the IRC02b ^i^^Ai?/ astronomical observing 
template (AOT) with the NP spectroscopy mode was used to obtain low resolution spectra [R 
~ 20) between 2.5 and 5 fim. By June 2008, 621 pointing observations were performed toward 
the LMC, and about 10 deg^ area was observed. The near-infrared (NIR) spectroscopic 
survey is continuing even after the exhaustion of the liquid helium in August 2008. 
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The spectral analysis was per formed using the sta ndard IDL package prepared for the 
reduction of AKARI IRC spectra (lOhyama et al.l 120071 ). Raw data were converted to dark- 
subtracted, linearity-corrected, and flat-field-corrected frames. Three pixels were integrated 
in the spatial direction to extract a point-source spectrum from the slitless spectroscopic 
image. Background levels of each source were estimated from the signal counts of the adja- 
cent regions on both sides of the target spectrum. T he wavelength calibration accuracy is 
estimated to be about ~0.01//m (lOhyama et al.l 120071 1 



3. The selection of YSOs 



We select infrared-brigh t objects from the point-source catalog of the Spitzer SAGE 
project ( Meixner et al. 20061 ) with the following selection criteria: (1)[3.6] - [4.5] > 0.3 and 
[5.8] - [8.0] > 0.6, and (2) [3.6] < 12 and [4.5] < 11.5, where [wavelength] represents the 
photometric va lue in magnitude at e ach wavelength in microns. Criterion (1) refers to the 
YSO model of IWhitney et al.l (120041 ). and criterion (2) comes from the detection limit of 
the AKARI IRC NP mode. This rough selection is applied to the sources located in the 
survey area of AKARI IRC, and about 300 sources are selected. These photometrically 
selected sources include not only massive YSOs, but also a large number of dusty evolved 
stars since their infrared spectral energy distributions (SEDs) are similar to each other. For 
the accurate selection of YSOs, we select the sources that show absorption features of the 
3.05 /um H2O ice and the 4.27 ^m CO2 ice stretching mode in their NIR spectra taken 
by the present spectroscopic survey. The presence of CO2 ice is strong evidence o f YSOs 
since the detection of CO? ice toward dusty evolved stars has not been reported (IBarlow 
I999I : ISylvester et al.lll999l ). Spectral overlapping with other sources located in a dispersion 
direction is a serious problem for slit-less spectroscopy, which makes it difficult to obtain 
reliable spectra. We check the overlapping contamination by visual inspection, and we only 
use the sources without such contamination in the following analysis. 

As a result, we spectroscopically confirmed seven massive YSOs in the LMC for the 
first time. The sources are listed in Table 1 with the observati on parameters. Six o f the 
seven sources are included in the recent YSO candidates catalog (IWhitney et al.ll2008l ). and 
one source (ST6) is a newly found YSO. The spectra of these sources are shown in Fig 1 
together with the results of spectral fitting (see §4 for details). The absorption features of 
H2O and CO2 ices are rather broadened due to the low spectral resolution of the AKARI 
IRC NP spectroscopy mode but clearly detected. This is the first clear detection of the 4.27 
/im stretching mode of CO2 ice toward extragalactic YSOs. In addition, unresolved emission 
of PAHs and the hydrogen recombination line Pf6 around 3.3 //m, the Bra line at 4.05 /xm 
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and blended absorption features of 4.62 /im XCN and 4.67 /im CO ices around 4.65 /im are 
detected toward several sources. However, it is difficult to evaluate the column densities of 
XCN and CO ices accurately with the present low spectral resolution data. 



4. Spectral Fitting 



We fit a polynomial of the second to fourth order to the continuum regions and divide 
the spectra by the fitted continuum (Fig 1). The wavelength re gions for the cont inuum are 
set to be 2.5-2.7 fim, 3.6-3.7 ^m, 4.0-4.15 /im, and 4.9-5.0 /im doibb et al.ll2004l ). 



Due to the low spectral resolution of the IRC NP spectroscopy mode, direct comparison 
of the observed spectra with the absorption profiles of laboratory ices is difficult. However, 
the equivalent width of absorption does not depend on the spectral resolution of the spec- 
trum. Therefore we used a curve-of-growth method to derive the column densities of the 
ices. A Gaussian profile with the fixed central wavelength is fitted to the absorption bands 
to derive the equivalent width (Fig 1). The vertical axis of the plot is shown in units of the 
normalized fiux because it is difficult to estimate the optical depth directly from the present 
low resolution data. We use the laboratory absor ption profiles of H2O and CO2 ices taken 
from the Leiden Molecular Astrophysics database (lEhrenfreund et al.lll996l ) to calculate the 
curve of growth. The profiles of pure H2O ice and mixture of H2O : CO2 (100:14) ice both at 
lOK are used for the calculation since th ese compositions are typical in the interstellar ices 
( iNummelin et al.ll200ll : iGibb et al.ll2004l ). The present spectrum cannot resolve the polar 
and apolar CO2 ice features, and the present analysis assumes the pol ar CO2 ice only. How- 



ever, contribution of the apolar ice is generally small toward YSOs (IGerakines et al.lll995 



I999I ) . The FWHM of the absorption profiles of H2O and CO2 ices change by approxirn ately 
10 % depending on the compositions of the ices (IGerakines et al.lll999l : lGibb et al.ll2004j ). We 
confirm that this range of the variation in the FWHM makes negligible effects on the derived 
column densities compared with the observational errors. We adop t the band strengths o f 
H2O and CO2 ices to be 2.0x10^^^ and 7.6x10"^''' cm molecule"^ (jCerakines et al. 1995 ). 
respectively. The derived column densities are listed in Table 1. 



Since this study applies the curve of growth method to the calculation of the ice column 
density for the first time, we check the validity of our method. We use a few Galactic 
YSO spectra take n by ISO SWS whose column density of H2O and CO2 ices are derived in 
Gibb et al.l (120041 ). These high-resolution spectra are converted to low-resolution spectra by 
convolving the slit function of the NP spectroscopy mode, and then we derive the column 
density for these converted spectra using the same method de scribed above. The comparison 
of the obtained column density with the value presented in iGibb et al.l (120041 ) shows that 
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the differences are within 10 % and 20 % for the H2O ice and the CO2 ice, respectively. The 
result plotted in Fig 2 includes this uncertainty in addition to the observational error. 



Results and Discussion 



The obtained column densities of H2O and CO2 ices are plotted in Fig 2. The error 
bars become larger for the larger column density due to the saturation effect of the curve- 
of-growth. A linear fit to the data points indicates that the CO2/H2O ice column density 
ratio in the LMC is 0.45 ± 0.17. The large uncertainty mainly comes from the errors in 



the curve-o: 
taken from 



-growth ari a lysis. For comparison, column densities of Galactic massive YSOs 



Gibb et all (j200J) and their CO2/H2O ice column density ratio of 0.17 ± 0.03 



( Gerakines et al. 19991 ) are also plotted in Fig 2. A si milar CO2/H2O rat io of 0.18 ± 0.04 is 



also observed toward a Galactic quiescent dark cloud (IWhittet et al.ll2007l ). while a relatively 
high CO2/H2O ratio of 0.32 ± 0.02 is observed toward Galactic low- and intermediate- mass 
YSOs, and some of them reach ~0.4 (jPontoppidan et al.ll2008l ). Although the uncertainty is 
large, it is clear from the present results that the CO2/H2O ice ratio in the LMC is higher 
than the typical ratios of the Galactic objects. Since the distribution range of the H2O ice 
column density in the LMC is comparable to that of the massive Galactic YSOs, it can be 
concluded that the abundance of the CO2 ice is higher in the LMC. The present results 
suggest that the different galactic environment of the LMC is responsible for the high CO2 
abundance. 

The formation mechanism of CO2 ice in circumstellar environments of YSOs is not un- 
derstood, however a number of scenarios have been proposed. Several laboratory experiments 
indicate that the CO2 ice is efficiently pro duced by UV photon irradiation to H2O-CO binary 



ice mixtures (e.g., IWatanabe et al.ll2007l ). The LMC has an order-of-ma gnitude stronger UV 



radiation field than our Galaxy due to its active massive star formation (llsrael fc de Graauw 



19861 ). which could lead to the higher CO2/ H2O ratio in the LMC. The high CO2/H2O ratio 



toward a YSO in the LMC is also reported in Ivan Loon et al.l (120051 ) . and they suggest that a 



different radiation environment in the LMC is one of the reasons for the high CO2 abundance. 
On the other hand, the model of diffusive surface chemistry sug gests that high abundance of 



CO 2 i ce can be produced at relatively high dust temperatures (IBergin. Neufeld. fc Melnick 



I999I : iRuffle fc HerbstI I2OOII ) . Several studies have reported that the dust temperature in 



Aeuirre et al. 


2003; 


Sakon et al. 


2006) 



servations of diffuse emission (e.g., 
high dust temperature may also have an effect on the high CO2 ice abundance in the LMC. 



It is difficult to separate the effect of the UV radiation field and the dust temperature on 
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the high abundance of CO2 ice in the LMC by our low-resoluti on NIR spectra. The 4.62 fxm 

XCN feature is known to be indicative of strong UV irradiation (IBernstein. Sandford. fc AUamandola 



2OOOI : ISpoon et al.ll2003l ). On the other hand, detailed profile analysis of the 3.05 /im H2O ice 



stretching mode and the 15.2 /xm CO2 ice bending mode should rev eal the temperature and 
compositions of the ices (lEhrenfreund et al.lll996l : lOberg et al.ll2007l ). Future observations of 
the XCN feature and the H2O and CO2 ice features with a sufficient wavelength resolution 
will be useful to investigate this problem. 



6. Summary 

We performed a near-infrared spectroscopic survey of the LMC with AKARI IRC. We 
spectroscopically confirmed seven massive YSOs that show absorption features of H2O and 
CO2 ices. This is the first detection of the 4.27 /im CO2 ice feature toward extragalactic 
YSOs. The derived ice column densities indicate that the abundance of CO2 ice is clearly 
higher in the LMC than our Galaxy. The relatively strong UV radiation field and/or high 
dust temperature in the LMC may be responsible for the observed high abundance of CO2 
ice. Our study shows the difference in the chemical composition around extragalactic YSOs, 
suggesting that extragalactic YSOs hold quite different environments from Galactic ones. 

AKARI is a JAXA project with the participation of ESA. We thank all the members of 
the AKARI proiect for their continuous help and support. We also thank Youichi Ohyama 
for great support in the analysis of spectroscopic data, and we thank the AKARI IRC team 
and the LMC team for their continuous support. This work is supported by a Grant-in-Aid 
for Scientific Research from the JSPS. 
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Table 1. OBSERVATION PARAMETERS AND COLUMN DENSITIES OF ICES 



Number AKARI ID Obs. ID Obs. Date Other Name RA[J2000] DEC[J2000] N(H20) N(C02)'' 

(IQl^ cm-2) (IQl'^ cm-2) 

STl J053931-701216 2211375 2007 Apr 12 05393117-7012166'' 05:39:31.15 -70:12:16.8 ^-^ti.l ^■'^-t'.l 

ST2 J052212-675832 2210229 2006 Nov 24 NGC 1936 05:22:12.56 -67:58:32.2 H-ltl'i ^■'^-I'.l 

ST3 J052546-661411 2213132 2007 Jun 22 05:25:46.69 -66:14:11.3 '^^.rtte IS-Sl"/ 

ST4 J051449-671221 2200073 2006 Jun 6 IRAS F05148-6715 05:14:49.41 -67:12:21.5 18.7l^J S-SlH 

ST5 J053054-683428 2213061 2007 Mar 13 IRAS 05311-6836 05:30:54.27 -68:34:28.2 31.714 5 12.4t^;8 

ST6 J053941-692916 2201137 2006 Oct 25 05394112-6929166^ 05:39:41.08 -69:29:16.8 Sg.lljslo 

ST7 J052351-680712 2210220 2006 Nov 29 IRAS05240-6809 05:23:51.15 -68:07:12.2 ... 55.3lin ? 



''2MASS ID 

''The source is in a cluster 
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Fig. 1.— AKARI IRC 2.5-5^m spectra of YSOs in tiie LMC. Left: Plots of logio flux(Jy) vs. 
wavelengtli(//m). Open circles represent the points that used for the continuum determination. 
Dashed lines represent derived continuum. Right: Plots of normalized flux (1 - F/Fc; F and Fc 
represent observed flux and continuum flux, respectively) vs. wavelength(//m). Shaded regions 
around 3.0/im and 4.3/im represent fitted Gaussian for absorption features of H2O and CO2 ices, 
and the areas correspond to the equivalent width. The positions of H2O, CO2, XCN, CO ice 
absorption bands, PAH emission bands, and hydrogen recombination lines are shown. 
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Fig. 2. — CO2 ice vs. H2O ice column density in a unit of 10 cm~ . Open squares with error 
bars represent the results of this study. Filled squares represent those of Galactic massive YSOs 
(jGibb et al.ll2004l ). Upper and lower solid lines represent CO2/H2O ~ 0.45 and 0.17, respectively. 
The sources ST6 and ST7 are not plotted due to their large errors. 



